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ABSTRACT 



Context. The central few 100 pc of galaxies often contain large amounts of molecular gas. The chemical and physical properties of 
these extragalactic star formation regions differ from those in galactic disks, but are poorly constrained. 

Aims. This study aims to develop a better knowledge of the spatial distribution and kinetic temperature of the dense neutral gas 

associated with the nuclear regions of three prototypical spiral galaxies, NGC 253, IC 342, and Maffei 2. 

Methods. VLA CnD and D configuration measurements have been made of three ammonia (NH 3 ) inversion transitions. 

Results. The (J, K) = (1, 1) and (2,2) transitions of NH3 were imaged toward IC 342 and Maffei 2. The (3,3) transition was imaged 

toward NGC 253. The entire flux obtained from single-antenna measurements is recovered for all three galaxies observed. Derived 

lower limits to the kinetic temperatures determined for the giant molecular clouds in the centers of these galaxies are between 25 

and 50 K. There is good agreement between the distributions of NH 3 and other H 2 tracers, such as rare CO isotopologues or HCN, 

suggesting that NH3 is representative of the distribution of dense gas. The "Western Peak" in IC 342 is seen in the (6,6) line but not in 

lower transitions, suggesting maser emission in the (6,6) transition. 



Key words. Galaxies: individual: NGC 253, IC342, Maffei 2 
lines: galaxies 



1. Introduction 
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A large fraction (~ 20% at z < 0.2, Brinchmann et al. I2004l > 
of all stars, and as a consequence many of the massive stars in 
our Universe, are being born in the central regions of starburst 
galaxies. The reservoirs for such starbursts are large concentra- 
tions of dense molecular gas, which in many cases are formed 
by the interaction of two or more galaxies, which triggers the 
infall of large quantities of material toward the central few hun- 
dred parsecs in these systems (Combes 2005). The typical sizes 
of the clouds formed in such environments are a few 10s of par- 
secs. Observations to-date suggest that the physical properties of 
the molecular clouds in these extragalactic environments are dif- 
ferent from Giant Molecular Clouds (GMCs) found in the disks 
of the Milky Way and other normal galaxies. The chemistry of 
such extragalactic circumnuclear clouds is comparable in rich- 
ness to the Orion or Sgr B2 star forming regions, as a line survey 
at 2 mm toward NGC 253 suggests (Martin et al. 2006b). The 
spatial densities within these clouds should be ^ 10 4 ctrT 3 , i.e. 
higher than in disk clouds, in order to counteract the tidal forces 
induced by a high stellar density and a supermassive central en- 
gine (Giisten 1989). Such densities have indeed been confirmed 
by a number of molecular multilevel studies (e.g. Mauersberger 
etal.lT995l WeiB et al. p2TjQ7l Costagliola et al. l20TTjl Aladro et al. 
1201 11 1. As in the circumnuclear region of our own Galaxy (Oka 
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et al. 120051) . these dense clouds may be embedded in a less dense 
(~ lOOctrT 3 ) and warm (~250K) molecular component. 

Ammonia (NH3) is the classical tracer of the kinetic temper- 
ature (Tidn) within the dense neutral interstellar medium (ISM; 
see Ho & Townes 1983). Unlike most other molecules, where 
the transition excitation depends both on the density of molec- 
ular hydrogen, n(H2), and kinetic temperature, 7\i n , the relative 
transition intensities of the inversion transitions of this symmet- 
ric top molecule depend mainly on T^n for a large range of 
densities. In addition, NH3 has numerous inversion transitions 
at centimeter wavelengths covering a large range of energies; it 
therefore probes a wide range of temperatures. While brightness 
temperatures of main isotopic CO transitions are indicators of 
Tkin only for nearby dark clouds, where beam filling factors are 
close to unity, NH3 multilevel studies can be used to determine 
kinetic temperatures even in cases where there is sub-structure 
smaller than the size of the beam. This is particularly useful for 
extragalactic sources. 

The first, and for a long time the only, extragalactic ammonia 
study was limited to IC 342, a nearby face-on galaxy with par- 
ticularly narrow spectral lines (Martin & Ho I1986t . Thanks to 
improved sensitivity and spectral baseline stability, several ex- 
tragalactic NH3 multilevel studies have been published during 
the past decade for the nearby galaxies NGC 25 3, M5 1, M82, 
an d Maffe i 2 (Henkel et al. 2000J Takano et al. 200QJ WeiB et 
al. 12001 at Mauersberger et al. [20031 Ott et al. 120051 Kinetic 



2 



M. Lebron et al.: Dense gas in nearby galaxies 



temperatures have been estimated toward the central regions of 
these galaxies on spatial scales of several hundred pc. Takano 
et al. d2005al l also presented Very Large Array (VLA) images 
of the NH 3 (1,1), (2,2) and (3 ,3) emission from NGC253, and 
Montero-Castano et al. (2006) presented a VLA (6,6) image of 
IC 342. NH3 absorption lines have been reported toward Arp 220 
by Takano et al. 2005b and, at intermediate redshift, toward the 
main gravitational lenses of the radio sources B0218+357 and 
PKS 1830-21 1 (Henkel et al. 120031120081 . 

In this paper, interferometric observations of the (J, K) = 
(1,1) and (2, 2) ammonia inversion transitions are presented for 
the central regions of the nearby galaxies IC 342 and Maffei 2, 
and of the (3, 3) transition toward NGC 253. The linear resolu- 
tions of these new NH3 measurements are ~ 50 pc. The primary 
goals are to determine the spatial distribution and to constrain 
the kinetic temperature of the NH3 emitting gas, and to relate 
these to the dominant physical and chemical parameters of the 
studied regions. 

2. Observations and data reduction 

The ( J, K) - (1,1) and (2, 2) inversion transitions of NH3, 
(rest f requencies: 23.6944955 GHz and 23.7226333 GHz, Lovas 
|2004| |. were observed toward IC342 and Maffei 2 on 22-27 
October 2001 using the D configuration of the Very Large Array 
(VLA) of the NRACfl The beam sizes are 3'.'8 x 3'.'4 at a PA. 
-85° for IC 342 and 3'.'6 x 3'.'2 at a PA. -89° for Maffei 2. We 
used a bandpass of 25 MHz centered at heliocentric velocities 
of +40 km s _1 for IC 342. The heliocentric velocity of Maffei 2 
is -40 kms -1 , but due to limitations in the LO tuning at the 
VLA for bandpasses of 25 MHz or more, we set the central ve- 
locity to -35 kms~'. The total bandwidth for both sources was 
divided into 15 spectral channels, each 1.5625 MHz wide (~19.7 
km s _1 at the observing frequencies), plus a continuum channel 
containing the central 75% of the total band. For the amplitude 
calibration, we used 0137+331, with an adopted flux density of 
1.05 Jy. The phase calibrators were 0304+683 with a flux den- 
sity of 0.43 Jy, and 0244+624 with a flux density of 1.6 Jy at 
1.3 cm. 

The NH 3 ( 3,3) transition (rest frequency: 23.8701292 GHz, 
Lovas [2004b was observed toward NGC 253 on 29-30 
September 2001 using the VLA DnC configuration. This config- 
uration was required because of the low elevation of this source. 
The resulting beam size was 279 x 2'.'2 at a PA. 58°. The obser- 
vations used a bandwidth of 50 MHz centered at a heliocentric 
velocity of 9 km s _1 and 16 spectral channels, each of them cov- 
ering ~ 40kms~ 1 . Since the systemic velocity of NGC 253 is at 
~ 230kms~ 1 , the transition is at the edge of the observed band. 
From a comparison with CO spectra (e.g. Mauersberger et al. 
2003) we estimate, however, that S= 5% of the flux is outside 
the observed spectral bandpass. Absolute flux density calibra- 
tion was obtained by observing 0137+331 (1.05 Jy). Phase and 
bandpass calibration were performed on 0120-270. 

The data were edited and calibrated following the standard 
VLA procedures and using the Astronomical Image Processing 
System (AIPS) developed by the NRAO. Observations from dif- 
ferent days were combined in the UV plane. We synthesized, 
cleaned, and naturally weighted maps of the observed regions. 
For each source, a continuum map was produced by averaging 
line-free channels. 

1 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 



In order to convert flux density per beam, S , to synthesized 
beam brightness temperature, 7b, we have used the relation 

S(Jy/beam)/r b = 8.18 10~ 7 6> a (asec)<9b(asec)(v(GHz)) 2 (1) 

(see Eq. 8.20 in Rohlfs & Wilson [T9"9"6j l, where ^ and 6> b ' indicate 
the Full Width at Half Power (FWHP) linewidths of the major 
and minor axes of the elliptical synthesized beam in arcseconds. 

3. Results 

3. 1 . Continuum emission 

Fig. lA.ll shows the A ~ 1. 3 cm continuum of IC 342, Maffei 2, and 
NGC 253. We estimate that the calibration accuracy of the con- 
tinuum is ~ 10%. The 23.7 GHz continuum emission of IC 342 
peaks toward a /2 ooo = 03 h 46 m 47 ! i85, (5 J2 ooo = 68°05'46"0 with a 
total flux density of 33mJy and a maximum of 7.1 mJy/beam. 
The emission is resolved with an apparent FWHP extent of 
972 x 576 (major x minor axis), which corresponds to a decon- 
volved source size of 874 x 570 at a PA. of +78°. This extent is 
compatible with the observations by Turner & Ho (1983), who 
separated the thermal from the non-thermal emission. The ther- 
mal 6 cm continuum emission determined by these authors is 
~40mJy for IC342. Considering this flux density and a spec- 
tral index a of -0.1 (S v oc v ff ), the expected thermal emission is 
~ 36mJy at 23.7 GHz, which is consistent with our measured 
value of 33 mJy. 

The 23.7 GHz continuum emission from Maffei 2 is centered 
at orj2ooo = 02 h 41 m 55 s .02, dnooo = 59°36'1772 with a total flux 
density of 40mJy and a maximum of 8.2mJy/beam. The emis- 
sion is resolved with an apparent FWHP extent of 976 x 570 (ma- 
jor x minor axis), which corresponds to a deconvolved source 
size of 97 1 x 374 at a PA. of ~ 0°. This extent is compatible with 
the higher resolution VLA 2 cm data by Turner & Ho (1994), 
who found that the extended continuum emission is dominated 
by non-thermal radiation. They assumed an a of -0.8 for the ex- 
tended continuum in Maffei 2. The non-thermal component is, 
however, spatially related to compact thermal sources. Taking 
the above mentioned spectral index, the expected flux density is 
31.8 mJy at 23.7 GHz, and is smaller than our value of 40 mJy. 
This suggests that a significant fraction of the 23.7 GHz contin- 
uum is due to thermal free-free emission. 

Toward NGC 253 we have detected 23.7 GHz contin- 
uum emission centered at crnooo = 00 h 47 m 33 s T6, 5j2ooo = 
-25° 17' 1775 with a total flux density of 468 mJy and a maxi- 
mum of 179mJy/beam. This flux density is compatible, within 
the estimated errors, with the data obtained by Geldzahler & 
Witzel £1983} and Ott et al. d20051 >. The emission is resolved 
with an apparent FWHP extent of 57 1 x278 (major x minor axis), 
which corresponds to a deconvolved source size of 472 x 177 at a 
PA. of 45°. The spatial extent of this emission is compatible with 
the higher resolution images obtained by Ulvestad & Antonucci 
(TT997b . 

3.2. The NH3 emission 

Our maps of the NH3 emission toward IC342, Maffei 2, and 
NGC 253 confirm that the single dish ammonia spectra previ- 
ously detected toward these galaxies (Mauersberger et al. 2003) 
are composed of the contributions of different giant molecular 
clouds (GMCs). In the analysis of the data, for each identified 
giant cloud, Gaussians were fitted to the ammonia spectra. Table 
Q]contains the fit parameters for each GMC for IC 342, Maffei 2, 
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and NGC 253, respectively. The average linewidths for IC 342 
are ~ 40 km s -1 while for Maffei 2 and NGC 253 they are larger 
by a factor up to four (see §3.2.31 ). This, and the low velocity 
resolution of ~20 km s _1 for IC342 and Maffei 2 and ~40 km 
s _1 for NGC 253 does not allow for a deconvolution of the line 
profiles into velocity sub-structures inside the GMCs or to de- 
convolve the satellite lines as it has been done for Galactic center 
clouds (McGary & Ho l2002l ). 

3.2.1. IC342 

Improvements in the performance of the VLA interferometer 
have resulted in a highly increased dynamic range for our im- 
ages of the ammonia emission from IC 342 as compared to Ho 
et al. (1990). Our channel maps of the ammonia (1,1) and (2,2) 
emission (Figure IA.2t reveal clear substructure not apparent in 
previous measurements. The ammonia emission is present be- 
tween 0.4 kms~' and 80 kms~' . 

Figure Q] shows the velocity integrated emission for both 
transitions. The ammonia emission is concentrated in four main 
giant clouds. The peak position of each cloud is given in Table 
Q] while the corresponding spectrum is shown in Figure Q] Each 
spectrum was fitted with a single Gaussian in order to obtain the 
central velocity and the apparent linewidth. 

In the regions selected around some of the GMCs (see 
Figure the overall velocity integrated fluxes of the (1,1) and 
(2,2) transitions are 0.54 and 0.30 Jy kms~'. This is significantly 
smaller than the total flux in these lines contained in our maps 
(1.43 Jy kms 1 for the (1,1) line and 0.89 Jykms 1 for the (2,2) 
line) and indicates that the regions selected represent only about 
35% of the total NH3 emission. Ammonia is therefore not re- 
stricted to the main identifiable giant molecular clouds but is also 
detected from the intercloud gas. Within the calibration uncer- 
tainties and accounting for potential differences due to different 
beam sizes, the measured total flux is similar to the single dish 
fluxes from Mauersberger et al. (2003) which are 1.1 Jy kms~' 
for the (1,1) line and 1.24Jy kms 1 for the (2,2) line. Our inter- 
ferometric data are therefore not missing significant amounts of 
single dish flux from the ammonia emission. 

3.2.2. Maffei 2 

The ammonia (1,1) and (2,2) spectral channel measurements of 
Maffei 2 are shown in Figure lA3l revealing velocity structure 
which forms four main giant clouds along a ridge extending from 
the north-east to the south-west (Figure|2|i. There is emission to- 
ward all 4 peaks. In addition, in peak B, which is close to the 
continuum emission peak, the NH3 (2,2) is seen in absorption 
against the 1.3 cm continuum in the velocity range between -20 
and 100 kms" 1 . Maffei 2 is located close to the Galactic plane 
and the line of sight material causes a visual extinction of 5.6 
mag (Fingerhut et al. 120071 ). The NH3 absorption we observe 
is most likely from gas within Maffei 2 and not from interven- 
ing gas in our own Milky Way since a) the observed lines are 
much wider than a few kms -1 and b) the velocities expected 
from gas in this quadrant of the Milky Way would be at neg- 
ative velocities. In Table Q] we give the results of Gaussian fits 
toward the peaks A-D. The fluxes of the (1,1) and (2,2) emission 
lines contained in the regions A — D are 1.3 and 0.8 Jykms 1 , 
while the total line flux from these transitions contained in our 
map is 1.8 and 0.8 Jy kms -1 . For the (1,1) line this indicates that 
the emission from regions A — D represents 70% of the total flux 
in the mapped area. Such a percentage may also be representa- 



tive for the NH3 (2,2) emission regions, but a direct compari- 
son is not straight forward due to the mixture of emission and 
absorption features in this line, which are difficult to disentan- 
gle. The single dish fluxes from the 40" Effelsberg beam are 1 .2 
and 1.1 Jykms 1 (Henkel et al. I20001 ). For the (1,1) line this is 
smaller than the total flux in our map, which is expected since 
the NH3 emitting region in Maffei 2 is not small compared to the 
Effelsberg beam. In particular, peak D is right at the half power 
level of the Effelsberg beam. 

3.2.3. NGC 253 

The left frame of Fig. [3] shows the velocity-integrated emission 
of the NH 3 (3,3) transition. With an overall extent of 40" x 5" 
along a PA. of ~ 45°, the NH3 emission is much more extended 
than the 1.3 cm continuum, which is located near the "center of 
gravity" of the NH3 emission. The upper right of Fig. [3] shows 
the integrated spectrum of the NHb(3,3) emission. The line shape 
and the flux of 8.6 Jy km s _I is only slightly smaller than the total 
flux detected in our maps (11.8 Jykms~'), which is similar to 
the single dish spectrum observed by Mauersberger et al (2003) 
with a 40" beam, which contains a flux of 10.9 Jykms -1 . The 
flux measured by Ott et al. (2005) with the ATCA interferometer 
(synthesized beam sizes 19" x 5" and 30" x5") is 8.1 Jykms -1 . 
Takano et al. (2005a) report a significantly higher value that is 
30% higher than ours, namely 15.3 Jykms -1 . This may be to a 
very small part due to our incomplete spectral coverage of the 
NH3 (3,3) line, but mainly reflects the calibration uncertainties 
at such a low elevation. We will assume in the following that, as 
for Maffei 2, most of the flux density has been recovered by our 
measurements. 

The NH3 emission arises from six giant clouds, that we la- 
bel from A to F with increasing right ascension. The lower-right 
frame of Fig. [3] shows the NHb(3,3) spectra observed toward 
each of the six giant clouds (indicated by boxes in the contour 
map). Each spectrum was fitted with a single Gaussian (see Table 
[TJ. Typical linewidths of the individual clouds are of order 70 — 
160 km s _1 . This is larger than the velocity separation of the hy- 
perfine (HF) components of the NH.3(3,3) transition. Under con- 
ditions of Local Thermodynamic Equilibrium (LTE) and opti- 
cally thin line emission, the outer hyperfine components have an 
intensity 2.6% of that of the main component. Therefore, the fit- 
ted linewidths should well correspond to the intrinsic linewidths 
of the clouds. Channel maps of the ammonia (3,3) emission are 
shown in Figure IA~4l The ammonia emission is present between 
144km s _1 and 341kms~ 1 . The latter is the upper limit to the 
velocity extent of the emission dictated by observational band- 
width limitations (see §|2](. 



4. NH , column density and rotational temperature 

The total column density of the two states of an ammonia (J, K) 
inversion doublet is given by 

N(J, K) = 1.65 lO^cnrV 1 [j(J + l)/K 2 ] AvT ex T tot , (2) 

where v is the inversion frequency in GHz, Av is the FWHP 
linewidth in kms -1 , T ex is the excitation temperature in K, and 
T tot is the sum of the peak optical depths of the five groups 
of hyperfine (hf) components comprising each inversion tran- 
sition. For the (J,K)=(l,l) transition, T tot =2.000r m (m: main 
group of hf components); and for the (2,2) and (3,3) transitions 
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Table 1. Gaussian fit parameters, NH3 column densities, and rotational temperatures 



Position 














a(J2000) (5(J2000) 


NH 3 


S 


FWHM 


v hel 


N(J, K) 


r rot (l, 1;2,2) 


(h m s) (• ' ") 


(J,K) 


(mJy) 


(km s" 1 ) 


(km s- 1 ) 


10 13 cnT 2 


K 


IC342" 














A 03 46 48.47 +68 05 43.2 


(1,1) 


1.9 


39.0 


23.6 


13.3 


35.8 




(2,2) 


1.3 


40.5 


21.5 


7.1 




C 03 46 49.10 +68 05 50.9 


(1,1) 


5.6 


44.0 


50.3 


16.1 


31.4 




(2,2) 


3.7 


40.2 


46.9 


7.3 




D 03 46 49.60 +68 05 58.6 


(1,1) 


3.1 


45.1 


52.2 


5.1 


24.7 




(2,2) 


1.4 


42.2 


50.4 


1.6 




E 03 46 47.47 +68 05 41.8 


(1,1) 


2.2 


37.4 


13.7 


16.5 


25.0 




(2,2) 


0.9 


39.4 


13.2 


5.3 




Total Map Flux 


(1,1) 

(2,2) 




1.43* 
0.89* 






29.3 


Maffei2 














A 02 41 55.66 +59 36 26.0 


(1,1) 


2.0 


65.6 


-82.6 


15.5 


29.0 




(2,2) 


1.3 


54.3 


-82.8 


6.25 




B' 02 41 55.20 +59 36 20.0 


1,1) 


4.8 


74.2 


-84.9 


17.0 


47.8 




(2,2) 


3.8 


88.3 


-85.9 


12.0 




C 02 41 54.89 +59 36 10.2 


1,1) 


5.6 


66.1 


0.15 


15.6 


23.9 




(2,2) 


4.0 


37.1 


-1.9 


4.7 




D 02 41 53.94 +59 36 59.3 


1,1) 


7.3 


57.2 


23.5 


9.5 


33.0 




(2,2) 


5.4 


49.6 


18.1 


4.56 




Total Map Flux 


(1,1) 
(2,2) 




1.83* 
0.78* 






32.8 
















A 00 47 32.05 -25 17 26.5 


(3,3) 


22.2 


85.0 


313.7 


74 




B 00 47 32.33 -25 17 20.0 


(3,3) 












C 00 47 32.81 -25 17 21.1 


(3,3) 


26.5 


62.2 


290.7 


122 




D 00 47 33.25 -25 17 16.6 


(3,3) 


17.8 


161.6 


210.9 


170 




E 00 47 33.66 -25 17 12.8 


(3,3) 


14.7 


79.1 


187.4 


123 




F 00 47 34.06 -25 17 11.2 


(3,3) 


14.8 


69.6 


214.9 


62 




Total Map Flux 


(3,3) 




11.8* 









" The nomenclature follows the one presented by Downes etal.t 1992 their Fig. 1). 
Component B is not seen in ammonia. 



* / Sdv(Jy Kkms- 1 ) 
c There is ~1 mJy of absorption in the NH 3 (2,2) line toward source B of Maffei 2 in the 
velocity range between -20 and 100 km s- 1 . Since this absorption is at the edge 
of our band we did not fit a Gaussian to it. 



Tt t=1.256T m , and 1.119r m , respectively. For higher-lying inver- 
sion transitions r tot - r m . Assuming that the background con- 
tinuum can be neglected, the brightness temperature of the main 
group of hyperfine components 7b of an ammonia transition re- 
lates to its excitation temperature as 7\, = T ex (l - exp(-T m )), 
i.e. in the optically thin limit 7b ~ T ex T m . Using the transition 
parameter that can be directly observed, namely the main-beam 
brightness temperature T m \,, instead of 7b results in beam aver- 
aged (instead of source averaged) column densities. 

The transition optical depths are unknown and cannot be eas- 
ily determined since Doppler broadening masks individual hy- 
perfine satellite emission. Takano et al. (2005b) claim low opti- 
cal depths from line-to-continuum ratios of absorption lines to- 
ward Arp 220. The latter is, however, also affected by the un- 
known source covering factor of the molecular gas. From ob- 
servations of ammonia emission toward molecular clouds to- 
ward the Galactic Center region, where linewidths are of or- 
der 1 km s~ 1 , Hiittemeister et al. ( 1 993 ) derive optical depths of 
about two for the main group of hyperfine satellites of the (1,1) 
and (2,2) transitions. They also find higher optical depths toward 
the stronger peaks of molecular emission. The Galactic Center 
data presented by Hiittemeister et al. Q993 ) were obtained with 
a linear resolution of 1.7 to 3.2 pc, while the VLA beam size to- 
ward the galaxies observed corresponds to about 50 pc (see ©. 



It is plausible that the source-averaged opacities of the observed 
NH3 transitions toward our galaxy sample are smaller than the 
value of ~ 2 determined toward Galactic Center clouds, which 
would suggest that T ex r tol ~ T mb . 

The synthesized beam averaged T m \, values are associated 
with the spectral flux density, S v , by the relation in EquationQ] 
With these assumptions, beam averaged column densities can be 
derived using 



N(J,K) = /dO^cm^XT'SjyAv-tkms- 1 ), 



(3) 



where k = 6.04, 2.81, 2.21, and 1.47 for the (J,K)=(1, 1), (2,2) 
and (3,3) and (6,6) transitions discussed in this paper, respec- 
tively. 

The rotational temperature between any (J, J) and (J', J') 
levels can then be determined using 



T mt = AE(K)/ln 



(2J+ l)g op N(J',J') 



(2J' + \)g' N{J,J) 



(4) 



Here g op is 2 for ortho ammonia (K = 3, 6), and 1 for para am- 
monia (K = 1,2). 
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Fig. 1. Shown as contours, emission integrated over the velocity range 0.45 to 59.8 km s~' for the NH3(1,1) transition (left) and 
NH3(2,2) transition (central) toward IC342. The contour levels are -3, 3, 6, 9 and 12 times 0.7 mJy km s beam -1 for the (1,1) 
map and -3, 3, 6, 9 and 12 times 0.5 mJy km s -1 beam -1 for the (2,2) map. Overlaid in gray scale for each panel are the radio 
continuum maps at the respective frequencies (~23.7 GHz) with the upper bars providing the flux density scale in mJy beam 1 . An 
asterisk indicates the location of the (western) water vapor maser detected by Tarchi et al. (2002 ), while the boxes indicate the areas 
that were taken for integrating the emission from each ammonia cloud identified from these data. The spectra to the right show the 
integrated NH3 emission (upper panel) and the spectra toward the four regions depicted in the contour maps (lower panel). 



5. Discussion 

5.1. IC342 

5.1 .1 . Distribution and excitation of the molecular gas 

Distance estimates toward the nearly face-on and highly ob- 
scured spiral galaxy IC 342 ran ge from 1.8 Mpc (McCall 1989) 
to 8 Mpc (Sandage & Tammann l 19741 ). As a consequence of this 
distance ambiguity the size and luminosity of this galaxy have 
been classified as similar to or much in excess of the Galactic 
values. Recently Tikhonov & Galazutdinova (2010) derived a 
distance of 3.9+0.1 Mpc using stellar photometry. To facilitate 
a comparison with previous studies, here we adopt the distance 
proposed by Saha et al. (2002), 3.3 Mpc. The giant molecular 
clouds in the central bar of IC 342 (e.g. Sakamoto et al. 1999 ) 
have a linear size which is similar to that of the Sgr A and Sgr B2 
molecular clouds in the center of our Galaxy. The inner 400 pc of 
IC 342 have the same far infrared luminosity as the inner 400 pc 
of our Galaxy and the 2fim luminosities indicate that the cen- 
tral regions contain similar masses of stars (e.g. Downes et al. 
11992) . Fortunately, IC 342 is viewed face-on, and its nuclear re- 
gion is therefore much less obscured than that of the Milky Way, 
although IC 342 is located behind the Galactic plane (b 11 ~ 10°). 

5.1 .2. Comparison of NH 3 with other tracers of dense gas in 
IC342 

According to Downes et al. d 19921 ) and Sakamoto et al. ( 1999), 
the highest concentrations of CO and HCN are present in a re- 
gion of diameter ~40" near the nucleus of IC 342. Downes et 
al. ([1992) discussed the presence of at least five molecular com- 
plexes, the most interesting of which is "Cloud B" (following 
the notation of Downes et al.). Only this region is associated with 



powerful H n regions, while the other complexes appear to be too 
hot and too turbulent to form stars with similar efficiency. The 
nucleus and Cloud B are only ~3" apart. The ammonia emis- 
sion is in the form of four unresolved or poorly resolved maxima 
along an elongated curved structure of 40" (350 pc) length. Four 
of these maxima can be identified with the HCN maxima A, C, 
D, and E observed by Downes et al. ( 1992). The ammonia ridge 
continues southward of HCN peak E without resolving single 
GMCs. No NH3 emission can be seen toward "Cloud B". 

When we compare our NH3 (1,1) and (2,2) emission with 
the various interferometric molecular emission maps by Meier 
& Turner J200T1 120051) and Usero et al. d2006i ). we find the best 
correspondence (including the extension south of cloud E) with 
the distributions of C 18 (1 - 0), CH 3 OH (methanol) (2 k - l k ) 
and N2H + emission. We also find excellent correspondence with 
images of H 13 CO + (1-0) and SiO (2 - 1), although in these 
Plateau de Bure maps the extensions south of Cloud E are not 
visible as they are outside the primary beam. There is also cor- 
respondence between HCN peaks A, C, D, and E and HNCO 
(4o4 - 3o3), although the extension south of Cloud E is not seen 
in H NCO. With HNC (1-0), as with HCN (1-0) (Downes 
et al. 119921 ), there is good coincidence of peaks A, C, D, and 
E; however, Cloud B is seen in HNC and HCN but (as already 
mentioned) not in NH3. The poorest coincidence with our NH3 
images is with the distributions of CCH 1 -0; 3/2- 1/2 and C 34 S 
2-1, both of which peak between Cloud A and the continuum 
peak. 

Of the more complex molecules that coincide the best, 
namely NH3, CH3OH and HNCO, all have in common that they 
are thought to be liberated from grains by shocks. In order to 
be released from grain mantles without being destroyed, slow 
(vshock = 10 - 15km s -1 ) shocks are required (see the discus- 
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Fig. 2. Images of the velocity-integrated NH3( 1,1) and NH3(2,2) lines from Maffei2. The emission was integrated from -135 km s _1 
to 62 km s _1 . The gray scale images are the continuum maps at the respective frequencies (~23.7 GHz) with the upper bars providing 
the flux density scale in mJy beam -1 . The contour levels are -3, 3, 4, 5, and 6 times 1.3 mJy km s _1 beam -1 for the (1,1) map and 
-3, 3, 4, 5, and 6 times 1.8 mJy km s _1 beam -1 for the (2,2) map. The spectra to the right show the integrated NH3 emission (upper 
panel) and the spectra toward the four regions depicted in the contour maps (lower panel). 



sion in Martin et al. I2006al . Ammonia, methanol and HNCO 
are easily destroyed by dissociating radiation (Hartquist et al. 
[19951 Le Teuff et al. EUOOl l. SiO, on the other hand would need 
fast shocks (v s h oc k ^ 15 - 29km s -1 ) to be released from grain 
cores (see the discussion in Usero et al. 120061 1 . Alternatively it 
has been proposed that SiO can be generated in a high tempera- 
ture environment in the gas phase (Ziurys et al. 1989). The good 
coincidence of SiO with the slow shock tracers methanol, am- 
monia and HNCO would support this latter interpretation. For 
all the mentioned tracer molecules, our peaks A, C, D and E can 
clearly be identified, and they all seem to avoid the continuum 
peak indicating the absence (or effective shielding) of dissociat- 
ing radiation. 

In contrast, C 34 S, CCH (and also to a certain degree HCN 
and HNC) are strongly enhanced toward peak B. In the case of 
CS, this can be explained by a chemical enhancement of the CS 
abundance of up to > 10~ 6 due to the presence of a small amount 
of ionized sulfur (Sternberg & Dalgarno ll995l l. Also the CCH 
abundance is predicted to be enhanced in photodissociation re- 
gions (see the discussion in Meier & Turner 2005), which would 
explain its existence toward the only peak associated with H 11 
regions. 

5.1.3. Kinetic temperature in the IC342 clouds 

In Table[T]we list the derived rotational temperatures between 
the (1,1) and (2,2) rotational levels for the four peaks detected in 
these transitions; values range between 25 K (Clouds D and E) 
and 36 K (Cloud A). Since the rotational temperatures between 
the (1,1) and (2,2) levels represent only a lower limit to the ki- 
netic temperatures it is interesting to compare these with those 
derived using higher excitation transitions. 

The VLA map of the ( J, K) = (6, 6) (E„ = 411 K ) tran- 
sition toward IC 342 by Montero-Castano et al. d20061 > allows 



for a direct comparison with our images of the lower excita- 
tion NH3, and thus to study the excitation conditions in the 
various giant molecular clouds of IC 342. Comparing the max- 
imum of the spectral flux density, S, in their (6,6) VLA map 
with Effelsberg data by Mauersberger et al. (2003), Montero- 
Castano et al. (2006) conclude that no flux is missing in their 
NH3 (6,6) image, i.e. the three peaks represent indeed 100% of 
the (6,6) emission. The situation is different if one compares the 
flux density integrated over the line profiles (i.e. J S dv). The to- 
tal flux in the VLA NH3 (6,6) image is 242 mJykms -1 , while 
the Effelsberg spectrum presented by Mauersberger et al. (2003 ) 
contains a flux of 330 mJy km s -1 , i.e. 27% of the NH3 (6,6) line 
flux is missing in the VLA (6,6) map, most of which is probably 
in the line wings. 

In Fig.|4]we show a version of our NH3 (2,2) map smoothed 
to the same angular resolution as the (6,6) ammonia map of 
Montero-Castano (2006). In the NH3 (6,6) map, which has a 
resolution of 7'.' 8 x 5'.'0, three peaks can be identified: The main 
NH3 (6,6) peak can be identified with our peak C within an ac- 
curacy of 1". The second NH3 (6,6) peak appears to be north 
of our peaks A and E by 3", south of the HCN peak B, and 
coincides with the 1 .3 cm continuum. It is remarkable that this 
second peak does not coincide well with most other molecular 
line distributions presented in Meier & Turner (2005), except the 
HNC ( 1 -0) and the HC 3 N (10-9) transitions. The third peak, the 
"Western Peak" observed in NH3 is the most mysterious since it 
does not coincide with any published molecular distribution, in- 
cluding our NH3 (1,1) and (2,2) images, nor with any prominent 
continuum features. A notable exception is a 1.3 cm water maser 
discovered by Tarchi et al. Q.QQ2 ; see 95. 1 .5ft . whose positional 
uncertainty is 5". The NH3 (2,2) peak D, on the other hand, has 
no counterpart in NH3 (6,6). 

In Fig.[4j3 we show Boltzmann plots toward the other NH3 
(6,6) peaks. For this comparison the (1,1) and (2,2) images were 
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Fig. 3. Left frame: Contours show the velocity-integrated NH3(3,3) transition emission from NGC 253 integrated from 144 km s -1 
to 341 kms -1 . Overlain is the gray-scale image of the continuum at 23.7 GHz with the upper bar providing the flux density scale 
in mJybeairT 1 . The contour levels are -3, 3, 6, 12, 18, 24 times 1.2 mJy km s~' beam -1 . Right upper frame: Integrated NH3 (3,3) 
spectrum. Right lower frame: Ammonia (3,3) transition spectra integrated over the six clouds of NGC 253 identified by the boxes in 
the velocity integrated NH3 (3,3) emission (left panel). 



smoothed to the same resolution as the (6,6) image. This has 
the consequence that there is also some weak (1,1) and (2,2) 
emission from the continuum peak. We can observe the same 
behavior as was noted in the various transitions measured with a 
single dish toward IC 342 by Mauersberger et al d20031 >. While 
the rotational temperatures between the (1,1) and (2,2) transi- 
tions are between 23 K (continuum peak) and 27 K (peak C), 
the (6,6)/(2,2) rotational temperatures are 680 K toward peak C 
and 920 K toward the continuum peak. These values are higher 
than the single dish rotational temperature of 440 K fitted by 
Mauersberger ( [20031) to the NH 3 (5,5), (6,6) and (9,9) data. In a 
warm environment, radiative transfer simulations predict higher 
T rot values for the higher-excitation metastable energy levels, 
gradually approaching Tr, even if the gas is characterized by 
a single kinetic temperature (e.g. Walmsley & Ungerechts 1983, 
Danby et al.[1988 Flower et al.[l995). Alternatively, there may 
be several gas components with different temperatures within 
our beam. 

It was noted in 33.2. II that within the selected regions our 
NH3 (1,1) and (2,2) transitions recover less than 50% and 75 % of 
the (1,1) and (2,2) single dish fluxes, respectively. This will not 
significantly alter the resulting rotational and/or kinetic tempera- 
tures derived from these measurements. The different molecular 
complexes we have discovered are not much larger than the syn- 



thesized beam width of the VLA. There may, though, be an ex- 
tended emission component within which the measured compo- 
nents are embedded. Furthermore, since the percentage of miss- 
ing flux is larger for the (2,2) than for the (1,1) transition, the 
rotational temperatures given in TableQ]are lower limits. 



5.1 .4. The nature of the "Peak C" and the "Continuum peak" 

The NH3 emitting gas in the central region of IC 342 should orig- 
inate from highly shielded regions since NH3 has a low dissoci- 
ation energy and can be easily destroyed by a strong energetic 
radiation field (Saha et al. |1983l ). The temperatures observed are 
comparable to the kinetic temperatures of <: 600 K found in ab- 
sorption toward the Hn regions in SgrB2 (Hiittemeister et al. 
|1995t Wilson et al. 2006). The hot gas in SgrB2 is of compara- 
tively low density (< 10 3 cm -3 ). Hiittemeister et al. ( 1995) con- 
clude that the gas toward Sgr B2 is probably heated by C-type 
shocks arguing that this kind of shock is able to liberate large 
quantities of NH3 from dust mantles and heat them to high tem- 
peratures without destroying them. Like in the case of Sgr B2 the 
CH3OH and SiO observed in IC 342 might be associated with the 
hot gas component, at least for peak C and the continuum peak. 
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Fig. 4. Upper frame: Contours showing our IC 342 NH3 (J, K) — 
(2, 2) image smoothed to the same resolution as the Montero- 
Castano et al. (2006) NH3(6,6) emission image. The gray scale 
shows the 1.3 mm continuum, and the asterisks denote the three 
(6,6) ammonia peaks. The H2O maser detected by Tarchi et al. 
(2002 ) is indicated as a cross, while the box indicates the area 
depicted in Figure [T] Lower frame: Boltzmann plot for IC342 
peakC and the "continuum" peak using the measured (1,1), (2,2) 
and (6,6) transitions of NH3. 



5.1 .5. The nature of the "Western Peak" 

Toward the western peak, no definite (6,6)/(2,2) temperature can 
be derived due to our non detections of the lower excitation tran- 
sitions of NH3. A similar case where there is a strong presence 
of NH3 (6,6) emission with a corresponding lack of emission in 
the lower excitation transitions of NH3 has also been observed 
toward our Galactic center (Herrnstein & Ho [2002), although 
on a much smaller linear scale. This lack of NH3 (1,1) or (2,2) 
emission toward an NH3 (6,6) emission peak has been explained 
by a geometric configuration whereby the lower-excitation tran- 
sitions experience foreground absorption by cool material along 
the line of sight. However, although this scenario might explain 
the absence of lower excitation NH3 transitions, it would be dif- 
ficult to explain why one does not see any other molecular emis- 
sion toward this position. Alternative mechanisms that can ex- 
plain the presence of (6,6) emission but not lower lying transi- 
tions include: 



- An overabundance of ortho ammonia (i.e. with K 

0, 3.6,...) by a factor of > 2 due to its formation in a very 
cold (i.e. << 23 K) environment. Unfortunately, this config- 
uration has never been observed in any source in the Galaxy. 
The hypothesis could be tested by observing the NH3(3,3) 
transitions toward the NH3 (6,6) peak. 

- A more viable hypothesis is that we are observing maser 
emission in the (6,6) transition but not in the (1,1) and (2,2) 
transitions. Such a maser as observed toward the Galactic 
star formation region NGC 6334 1 (Beuther et al. 2007) could 
be detectable even if the column density is comparatively 
low as the detection of 15 NH3 ( J, K)={3,3) maser emission 
toward a Galactic hot core suggests (Mauersberger et al. 
1986). Ammonia maser activity has also been proposed to- 
ward NGC 253: Ott et al. d2003l) reported NH 3 (3,3) in emis- 
sion toward a position where other, higher and lower lying 
transitions were seen in absorption. In order to prove or dis- 
card the hypothesis of NH3 (6,6) maser activity toward the 
western peak, one has to perform higher resolution mapping 
and/or to study other high excitation transitions of ammonia. 

It is interesting that the western peak coincides, within the 
accuracy of the measurements, with narrow (Av ~ 0.5kms~') 
water maser emission detected by Tarchi et al. (120021 ). This water 
maser was found to be time variable on scales of weeks. Its inte- 
grated luminosity of Lh,o = 5 1O~ 3 L is comparable to that of 
the SgrB2 star forming region near the center of the Milky Way. 
While the H2O flux from Sgr B2 is spread over many velocity 
components, only a narrow flaring feature was observed toward 
IC 342. This is likely a consequence of its much larger distance, 
only revealing an occasional flare while sensitivity limitations do 
not allow us to detect the more "quiescent" components. Tarchi 
et al. (120021 ) found that toward the location of the IC 342 water 
maser (and, hence, toward the "western peak") there is a chain 
of optically identified sources that appear to be H 11 regions. 

5.2. Maffei2 

5.2.1 . Distribution and excitation of the molecular gas 

The barred, highly inclined (i = 61°) spiral galaxy Maffei2 
at a distance of ~ 3.3 Mpc (see the discussion in Meier et al. 
2008) is one of the nearest starburst galaxies, and hence, can 
be studied with high spatial resolution. A moderately strong 
nuclear starburst has been observed in the infrared (Rickard & 
Harvey 1983 ) and in radio continuum emission (e. g. Turner & 
Ho |1994| ). Its nuclear star formation rate is similar to those of 
M 83 and NGC 253 (Ho et al. I1990al ). but its low Galactic lat- 
itude (b u ~ -20 ) makes optical studies difficult (e. g. Buta & 
McCall [T9"83T ). 

Maffei 2 is a galaxy with strong molecular emission from its 
central region. CO emission reveals a ~ 220 pc long bar (Meier 
et al. 120081 ). which is distinct from the bar seen at large scale in 
this galaxy. From the location of dense gas and star formation 
near the intersections of X\ - X2 orbits, Meier et al. (2008 ) con- 
cluded that the starburst in Maffei 2 is triggered by dynamics. 

5.2.2. Comparison of NH 3 emission with other molecules 

Four molecular cloud complexes (designated in our Fig. [2] as 
NH3 A-D) can be seen in NH3 (1,1), two of which are situated 
north of the nucleus and the other two toward the south. The 
NH3 (2,2) emission is mainly seen toward the south of the nu- 
cleus. The large scale distribution of the NH3 (1,1) emission in 
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Maffei 2 follows that of the i sotopic CO and the HCN (7=1-0) 
emission (Meier et al. 120081 . 

If one goes into a detailed comparison of the various mea- 
sured molecular distributions, one can see important differences 
between NH3 and the other molecules. In 12 CO (1 - 0), and also 
in 13 CO (1 - 0), the emission is not only concentrated toward 
molecular cloud complexes but also a more widespread diffuse 
component can be seen. This might arise due to a higher opacity 
and a better signal to noise ratio in these transitions. The rarer 
isotopic substitution C 18 in its J = 1-0 and 2-1 transi- 
tions indicates, on the other hand, a similar contrast as NH3. The 
2-1 transitions are reliable tracers of the H2 column density 
for a large range of H2 densities (Wilson & Mauersberger 1990). 
In 13 CO 2-1, we can well distinguish the four NH3 clouds. 
Unlike NH3, CO also shows strong emission toward the contin- 
uum peak. The detailed distribution of C ls O 2-1 differs from 
the distribution of the NH3 and the 13 CO emission. This could 
be explained in terms of large variations of isotopic abundances 
of 18 or 13 C, contamination by the continuum, or just insuffi- 
cient signal-to-noise (D. Meyer, priv. comm.). To conclude, the 
similar distributions of large scale NH3 emission with other trac- 
ers of H2 suggests that we are indeed observing a representative 
fraction of the dense gas in NH3 and that the observed NH3 dis- 
tribution is not just due to excitation effects or due to varying 
chemical abundances. 

5.2.3. Kinetic temperatures 

In Table[T]we show rotational temperatures between the (1,1) 
and (2,2) levels for the four peaks detected in these transitions. 
Derived rotational temperatures range from 24 K (Clouds C) to 
48 K (Cloud B). Since the T m between the (1,1) and (2,2) lev- 
els represents only a lower limit to 7k, it is interesting to com- 
pare these measurements to rotational temperatures derived from 
higher-excitation NH3 transitions. Unfortunately such higher 
transitions are only available from single dish measurements 
(e.g. Henkel et al. [2000, Mauersberger et al. 120031 Takano et 
al. |2005a). Radial velocity information, can, however, be used 
for a comparison. From the (1,1) to (4,4) transitions, Henkel 
et al. (120001) suggested that the -80 km s _] component (Clouds 
A and B) should be warmer than the component at +6 km s -1 
(Clouds C and D). This is also what is indicated from our high 
resolution data. On the other hand the inclusion of (6,6) data 
by Mauersberger et al. (2003 ) shows no difference in excitation 
for the very hot gas. This indicates that although there is very 
warm gas (130K) at -80 and at +6 kms -1 , the gas toward the 
velocity component at +6 kms _I , i.e. clouds C and D has, on 
average, lower than average excitation, which could be due to an 
additional cooler component. 

5.3. NGC253 

5.3.1. The distribution of NH 3 and other molecules 

The Sculptor galaxy NGC 253 is a highly inclined (i ~ 78°) 
nearby (D ~ 3.5 Mpc, Recola et al. (2005 ) barred Sc spiral. It 
is one of the brightest sources of far infrared emission beyond 
the Magellanic Clouds and has been studied in detail at many 
wavelengths. NGC 253 hosts a nuclear starburst. Radio contin- 
uum measurements indicate that NGC 253 contains a large num- 
ber of potential centrally located supernova remnants and H 11 
regions (Ulvestad & Antonucci l 19971 Brunthaler et al. 2 0091 1 . 

A nuclear bar has been imaged in the CO, HCN and HCO + 
emission lines (Sakamoto et al. 120061 Knudsen et al. 2007 ), and 



many other molecules have been detected in NGC 253 (see e.g. 
the 2 mm line survey of Martin et al. [2006b ). Among these, five 
turned out to be particularly important to reveal the peculiar 
physical conditions of the starburst region in this galaxy: SO, 
SiO, HNCO, CH3O H and CH3CN (Henkel e t al. ITW71 [199111 
Nguyen-Q-Rieu et al. |199U Huttemeister et al. |1997l l. The abun- 
dances of these five molecules have been found to be particularly 
high in NGC 253. This indicates that in contrast to a "late star- 
burst" such as M82, where the chemistry is affected by pho- 
todissociation, the chemistry in NGC 253 is dominated by slow 
shocks which are able to release molecules from grains without 
destroying them (see also Garcfa-Burillo et al. 1200 II Usero et al. 
2004, 2006 ). Ammonia belongs to the same class of molecules, 
thought to be enriched in C-shock environments and which is 
easily destroyed by dissociating radiation due to its low dissoci- 
ation energy (Suto & Lee ll983l ). 

Takano et al. (J2005a) presented VLA maps of the (1,1), (2,2), 
and (3,3) transitions of NH3. The data presented in this paper 
only include the ammonia (3,3) transition but with a higher spa- 
tial resolution (2'.'9 x 2'.'2 instead of 4'.'5 x 2'.'5). The agreement 
between the two (3,3) maps is good. Comparing this map with 
a super-resolved NH3 map (resolution: 5" x 5") by Ott et al. 
(2005 ) we find excellent agreement with the locations of clouds 
A, B and C, while clouds D, E and F are shifted between 2 to 
5" with respect their peak positions. This discrepancy can prob- 
ably be explained by the deconvolution process used by Ott et 
al. (12005'). There is also good agreement between the features 
seen in our map and the 3'.'8 x 2'.'6 resolution maps of HCN and 
HCO + presented by Knudsen et al. (120071 ) 

The individual GMCs in NGC 253 are more pronounced in 
NH 3 than in HCN or HCO + . This may be due to the fact that the 
latter transitions being optically thick, and therefore the contrast 
between the clouds and the intercloud gas does not reflect the 
contrast in column density of HCN or HCO + . Another notable 
difference is that our Cloud B is about 3" north of the nearest 
HCN and HCO + feature. This HCN and HCO + feature, where 
no NH3 is measured, corresponds in location and velocity to one 
of the two "molecular superbubbles" identified by Sakamoto et 
al. (f2006 ), and is within the positional uncertainties identical to 
the superbubble already seen in the NH3 data of Ott et al. (2005 ). 
Also toward the other superbubble seen in the CO data presented 
by Sakamoto et al. (2006) and being located in the northeastern 
edge of the measured NH3 distribution (Fig. [3), weak HCN and 
HCO + emission can be seen at the correct velocity, but no NH3. 
At the opposite edge of the NH3 distribution, our Cloud A is 
clearly visible in HCN, but not in HCO + . 

The agreement of the NH3 distribution with that of 12 CO 
(2 - 1), 13 CO (2 - 1), 13 C O (3 - 2) and C lg O (2 - 1) measured 
by Sakamoto et al (2011) is excellent. Taking into account our 
sightly coarser resolution we can reproduce all the details in the 
distribution of those transitions. Since CO isotopes are thought 
to be among the most reliable tracers of the H2 column density, 
this supports our notion that NH3 traces the bulk of the gas in 
the central molecular zone of NGC 253. The reason may be that 
at the high temperatures in this region the abundance of NH3, 
which otherwise can vary by several orders of magnitude due 
to freeze out onto grains, is quite constant, since all NH3 is in 
the gas phase. The same may be true for the central molecular 
zone of the Milky Way, which has similarly high temperatures 
to a degree that even SiO can be observed throughout the region 
(Riquelme et al. l20T0l . 

To conclude, our NH3 distribution follows closely that of 
C 18 (2-1), which is thought to be an excellent indicator of 
molecular hydrogen. This suggests that the temperatures derived 
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from NH3 are representative of the bulk of the molecular gas in 
the central region of NGC 253. HCN and HCO + also follow the 
general H2 distribution. However, their emission is, with respect 
to NH3, enhanced toward regions where fast shocks and/or pho- 
todissociating radiation can be expected. 

5.3.2. Ammonia maser emission or unusual ortho/para 
ratios? 

Takano et al. (2002) noted that the (3,3) transition was stronger 
than expected toward NGC 253, and concluded that an anoma- 
lous ammonia ortho/para ratio was the cause. Also Ott et al. 
(2005 ) explained the unexpectedly high intensity of the (3,3) 
transition toward the NE peak in terms of an ortho/para ratio 
for ammonia of 2.5-3.5. On the other hand, Ott et al. (2005) 
also mention that toward the continuum position of NGC 253 
they measure the (3,3) transition in emission, whereas all other 
transitions are observed in absorption. This would suggest that 
maser emission in the (3,3) transition is the cause of the higher 
than anticipated intensities, which is amplifying the background 
continuum. 

In our measurements of the integrated (3,3) spectrum (see 
Fig. [3]) the high velocity SW component (V ~300 kms~') 
is brighter than the lower velocity NE component (V ~160 
kms~'). This is contrary to the lower resolution data presented 
by Ott et al ©05) and Mauersberger et al d20031 > where the 
170 kms" 1 component is equal in intensity or stronger than the 
300km s component. Note that consistently in all other am- 
monia transitions observed by Mauersberger et al. (2 0031 ) and 
Ott et al. (2005 ) the 300kms _1 component is the strongest. The 
flux density of the 170km s" 1 component is comparable in our 
measurements and those presented by Ott et al. The 300km s _1 
component, on the other hand, is 50% stronger in our measure- 
ments. From the channel maps in Fig.[3]and also from Takano et 
al. d2005ab . it appears as if this extra flux could come from cloud 
C (clouds SW1 in the notation of Takano et al. 2005a) which is 
close to the continuum peak. 

How can the different profiles be explained other than by cal- 
ibration uncertainties and errors? One possibility is time variable 
emission on a time scale of a few years. If the (3,3) transition is 
masing this would imply that the emission comes from a region 
small in angular extent near the center of NGC 253. Such maser 
emission would also imply that the intensity is very high within 
a very narrow transition linewidth. 

Inconsistent sensitivity to spatial structure could be an al- 
ternative explanation. The situation here is that the higher reso- 
lution VLA map results in a higher line flux than that derived 
from the Australia Telescope Compact Array (ATCA) image 
presented by Ott et al. (2005). If there is some large scale ab- 
sorption component on top of the emission from the GMCs, this 
might be "missed" by the VLA but not by the ATCA and may 
explain why the 300 km s component is stronger with the VLA 
than with ATCA. Apparently higher spatial and spectral resolu- 
tion observations not only of the (3,3) transition but also of other 
ortho or para transitions are necessary to settle this issue. 



6. Conclusions 

Imaging measurements of the NH3 (J, K) = (1, 1) and (2,2) tran- 
sitions toward IC 342 and Maffei2 and of the NH3 (3,3) transi- 
tion toward NGC 253 have been presented. Derived lower-limits 
to the kinetic temperatures determined for the giant molecular 
clouds in the centers of these galaxies are between 25 and 50 K. 



In general, there is good agreement between the distributions of 
NH3 and other H2 tracers, suggesting that NH3 is representative 
of the distribution of dense gas in these galaxies. 

For IC 342, a comparison to high resolution NH3 data indi- 
cates a molecular component with kinetic temperatures as high 
as 700 — 900 K. Furthermore, the "Western Peak" of IC342, 
also known to host an H2O maser, is observed in the NH3 (6,6) 
transition but not in lower-excitation NH3 transitions. This is 
suggestive of maser emission in the (6,6) transition. 

Comparing the (3,3) line profile from NGC 253 with 
those obtained from earlier lower resolution studies, the V 
~300kms 1 component appears to be enhanced. Explanations 
could involve a time variable (3,3) maser or a large scale absorp- 
tion component not seen by us. 
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Fig. A.l. Continuum maps at 23.7 GHz for IC 342 (left), Maffei 2 (middle), and NGC 253 (right). The contour levels are, -0.3, 0.3, 
0.6, 1.2, 2.4, 4.8 mJy beam" 1 for IC342, -0.6, 0.6, 1.2, 2.4, 4.8 mJy beam" 1 for Maffei 2, and -2.7, 2.7, 5.4, 10.8, 21.6, 43.2, 86.4 
mJy beam 1 for NGC 253. The beams are shown at the bottom left corner of each map. 
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Fig. A.2. Channel maps of the NH3(1,1) and NH3(2,2) transition emission from IC 342. The contour levels are -0.9, 0.9, 1.8, 2.7, 
and 3.6 mJy beam -1 and -0.75, 0.75, 1.5, 2.25, and 3 mJy beam -1 , respectively. The synthesized beams are shown in the lower left 
corner of the first panel for each transition. The gray scale images represent the continuum emission at 23.7 GHz. 
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Fig. A.3. Channel maps of the NH 3 (1,1) and NH 3 (2,2) emission from Maffei 2. The contour levels are -0.78, 0.78, 1.56, 2.34, and 
3.12 mJy beam -1 and -1.2, 1.2, 2.4, 3.6, and 4.8 mJy beam -1 , respectively. The gray scale images show the continuum emission at 
23.7 GHz. Beams are shown in the left upper panel for each map set. 
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Fig.A.4. Channel maps of the NH 3 (3,3) transition from NGC253. The contour levels are -1.8, 1.8, 3.6, 5.4, 7.2, and 9.0 
mJy beam -1 . The synthesized beam is shown in the left lower corner in the first panel. The gray scale images represent the continuum 
emission at 23.7 GHz. 



